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The Yakutsk BAS array has been designed for detecting 
the showers generated by the l0—10 eV primary cosmic 
rays and consists of numerous eleotron t muon, and Cerenkov 
light detectors arranged on a 20-km® area terrain (see Fig, 
4), The array is featured by the feasibility to detect the 
BAS-produeed Cerenkov light, henoe, as will be shown below* 
to find the mean energy of the primary particles generating 
an ensemble of EAS of given size. 

The Yakutsk array detectors are on the average spaced 
a relatively email distance apart (con. pared, for example, 
with the Havorah Park and University of Sydney arrays), the- 
reby permitting a comparatively high accuracy in determi- 
ning the basic parameters of the EAS detected. For instance, 
the Monte Carlo calculations allowing for the fluctuations 
of t e lateral distribution function of charged part ides 
have shown til that (see Table 1) the error in finding the 
parameter of individual SIS* which is the partiole 

density at a 600-m distance from the EAS axls*\ does not 
exceed 2 % if the shower axes fall within the effective de- 
tection area of the Yakutsk array for EAS of a given size. 
The calculations were oarried out for the various voltes of 
P gQQ and the various zenith angles of the EAS axes. The ef- 
fective detection area Is determined in terms of the requi- 
rement that the probability of EAS deteotion within the 

— tfT mm ~ - 

'The value of J* gg 0 can be determined most accurately 
with the Yakutsk array. Besides, fggo is known to be 
a good measure of the primary energy. 
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area should bo at least £ ■ 0,9, Irrespectively of the 
possible fluctuations in the form of the lateral distri- 
bu ion function. Fig. 2 shows the effective detection area 
for FA 3 with F 0 ^10^eV. The primary energy spectrum is 
inferred from the Yakutsk array data by finding, first, 
the BAS ^gQQ spectra for various intervals of the zenith 
angles of the detected BAS axes. The EAS axis direction 
and position end the EAS size ( f are found by the ma- 
ximum likelihood method described in d: tail in [l ] , 

Fig. 3 shows the BAS f goo s P ectra for the various ze- 
nith angles corresponding to the traversed atmospheric 
depths of 1046 g/cm 2 , 1133 g/cm 2 , and 1313 g/ cm 2 [l], The 
data of Fig, 3 may be used to find the experimental absor- 
ption path of EAS with a given value of j^goo* ^ 52 s/cm 2 
which, in turn, may be used to scale the observed spectra 
to the vertical direction $ * 0° (it should be noted that 
the same procedure is used to infer the .fg^Q spectrum 
from the Ilaverah Fork array dote) * 

Fig. 4 compares between the f g 0Q spectra inferred 
from the Yakutsk data within 3800® hours and from the Ha- 
verah Fork data. The I; aver ah Perk spectrum has been obtai- 
ned by scaling ?5Q0HP^ rom water— tank data to the 5- 

-om thick scintillator readings using he formula [2) • 

0 » (1 72-»0 25) f 1,06—0,03 

J 600 Sc J 600 HP 

on the basis of the 5* 600 HP 0 P Go ^ ram data from [3] and the 
J 7 goo spectrum data of the Yakutsk array from [l] . From 
Fig. 4 it Is seen that the f goose ®P 0c * 3Pa measured at Ya- 
kutsk and Havorali Park are in a good a reement with each 
other at goose ** ' ® n< * somewhat at higher 

y goose* ‘ ieceuse of ® limited accuracy in measuring the 
portiole flux densities with th* Yakutsk array, the effect 
of its geometry (for example, at the boundary of the array) 
the limited accuracy in finding, the axis position end di- 
rection and the partufeeler £ g 00 * the inaccuracy in calcu- 
lating the effective detection areas, and the detected EAS, 
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selection effect, the necessity arises to test the distort 
ting role of all these factors in terms of the Konte Carlo 
method, that is, to find the 5* g 0 0 0 P eo ^ rum distortion ‘fun- 
ction, She distortion function was calculated in [l] on tho 
assumption of the a priori power law £ gQOso spectrum Ftyj'-y 
allowing for the real errors in determining t e densities 
with the Yakutsk array end for the real aeleotion of the 
EAS to be studied* From Fig, 5 it is seen that the distor- 
tion function differs little (by not more than 10-2Q55)from 
unity at ^Vm 2 ) end ** IGd/m 2 ) , The corrections at 

y "'I and 10(1/m 2 ) arise from the effect of the EAS selec- 
tion system on measuring the particle flux densities with 
soincillators if the latter are included in the system. At 
f >20(i/m ), the corrections are even smaller. The EAS 
S 600 a P 4ct;ruffi may be used to infer thejenergy spectrum of 
primary cosmic rays from the experimental time-integrated 
arid differential Cerenkov 11 ht flux in an IAS ensemble 
with a given value of J* g 00 , As was first shown in [4] p the 
integral EAS Cerenkov light flux Q is directly related to 
the energy lost; by the s ower particles for ionisation 
and for the excitation of the medium atoms above the EAS 
observation level, namely, E^ ** k<3, where, according to the 
calculations [5], k depends little on the position of the 
EAS maximum. On the other hand, the data on the time-dif- 
ferential Cerenkov light flux at great distances from EAS 
axis make it possible to find the position of the maximum 
of an individual EAS [ 6 ] and, hence, to ^Judge about the ex- 
tent to which is close to the primary energy E ,• In the 
range of y g 0() » 1-IOffCI/m 2 ), considering the altitude of 
the maxima of the respective showers, we obtain that the 
factor k in the relation « kQ is 3»dxi0* eV/photon/eV to 
within several per cent. At the same time, because of the 
hi h location of EAS maximum for f qqq * 1-l00(i/m 2 ), it 
appeared that ® 0,8E Q ; in other words, the energy scat- 
tered by a shower above the observation level makes the ma— 
5 or contribution to the energy of the primary particle go- 
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aerating an RA3 ensemble with a given f 600 [l] ** , tmen 
finding the value of B^ (which la the major component of h Q ) 
it is absolutely neeesnary to allow for the EAS Cerenkov 
light absorption in the atmosphere due to the Rayleigh 
scattering and the scattering by aerosols ♦ The value of § 
measured directly may differ from the number of the gene- 
rated photons involved in the relation «* IdJ, 

Bearing in mind the above/ wet carried out the regular 
control of the atmospheric transparency above the Yakutsk 
array using the large optical detector described in [$0 • 

The occurrence frequency of the Cerenkov light bursts in 
the detector exceeding a certain threshold is contingent 
on the atmospheric transparency at the moment of measure- 
ment and on the spectral index of the burst intensity speo- 
trum E ( > q) •* (q * T) * where T is the atmospheric transparen- 
cy for the Wavelengths studied by the detector} q is the 
burst intensity; x is the index of the burst spectrum. 

Fig *6 shows the distribution of the burst occurrence 
frequency inferred from the 15-min observation intervals 
during several months of the Yakutsk array operations. The 

A 

detected burst threshold is 17 photon/eracV, The maximum 
value of K = 400 , If the variations of the Cerenkov 

light flux are assumed to be due only to the changes of 
the atmospheric transparency, the distribution shown in 
Fig, 6 may bo transformed into the distribution of S with 
the mean value T * 0*71T^ and ^(f )/§ * 0,19, where is 
the maximum transparency which corresponds to the Rayleigh 
scattering disregarding the aerosols. Considering the an- 
gular distribution of the axes of the EAS detected with 
the Yakutsk array, one may calculate [ 9] the value of 

*^In [l] it has been shown that, if the mean energy of EAS 
muons is taken to be E^ » 9 GeV according to the expe- 
rimei.tal data C7] and the lateral distribution of elec- 
trons near EAS axis is determined by the I3KG function 
with S » 1,15 C®]» the enery carried by electrons and 
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which proved to be Q.ii5, Since T « 0.71T,j, then T « 0.6. 

Thus, the true value of Q dif ers from its measured 
value by the factor 1/T ^*1,65. It should be noted that a 
a more rigorous allowance for toe atmospheric transparency 
involves also the usage of the ratio ^ ( T ) /T in making the 
scaling. However, the value of the ratio is suf ioiently 
small and, to within several per cent, does not affect the 
results. It should also be noted that t e large optical de- 
tector is a looal selection system which gives an idea of 
photon absorption by a layer of several km thickness above 
the observation level C9^ . The coincidence -mo d e functioning 
of the Samarkand array optical detectors spaced §20 meters 
apart has shown that Cio] the number of coincidences in the 
space d-epart detectors is in a good correlation with the 
number of coincide; cos in the local optical detectors at Sa- 
markand <K «* 0.9). This fact indicates that the EA8 Ceren- 
kov light absorption occures mainl. within the layer imme- 
diately above the observation level, rather than uniformely 
throughout the atmosphere. This conclusion coincides with 
the present-day concepts concerning the atmospheric aerosol 
layer of a * 2 -km thickness fill * 

Sc, considering the correction for the atmospheric 
transparency, we may obtain [z] thelfollowing experimental 
relation between the value of f qqq in KAO and the mean 
energy E of the primary particles responsible for an EAS 
ensemble with a given f qqq* 

B q * (5.0-1.4)3d0 l7 J , 600 0 * 96t °* 04 ( $ = 0°), 

The error in the nume< icol factor is mainly due to the er- 
ror (25, ) In the absolute calibration of the Cerenkov de- 
tectors. 

The statistical errors, which are of importance when 

muons below he observation level will be « *\% of R , 

The rest energy belongs to neutrinos and to nuclear 
splitting®. 
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deriving the above relation, can be seen in Fig,7. A given 
value of f g 0 Q la, generally, In correepondence to a dis- 
tribution of How ever, the width of the E 0 distribution 
is 00 small ( /D(E^)/E 0 0.2) that the relation between f 0 

and ^goQ is quite sufficient to use Ct© within several per 
cent) in making the transition from the £ g Q0 spectrum to 
the primary particle energy speotrum. 

Fig ,8 shows the experimental differential energy spec- 
tra of primary cosmic rays inferred from the Yakutsk array 
data Cl,2l, from the Haverah Pork data [j}, and from the 
Fly's Eye data Ol2l. Fig. 9 shows the experimental energy 
spectra obtained also in [l,d,12], From Figs 8 end 9 it is 
seen that at energies below IQ^eV the Yakutsk and Haverah 
Park data are in a better agreement with each other than 
with the Fly's Eye data. The deviations of the latter are 
probably due to the tentative oharaoter of the results of 
calculating the Fly's Eye geometric factors at low primary 
energies, as was noted in [l2l. The minor disagreement bet- 
ween the Yakutsk and Haverah Park data arose probably from 
the foot that the Hillas model used for the Haverah Park 
array to scale from fg 00 to E q is not sufficiently accu- 
rate in reflecting the true relation between f gQO B o* 
At E c > lO^eV the experimental data of different works show 
a significant spread. Such a spread may be caused by tri- 
vial reasons, for example, e certain inaccuracy in the ab- 
solute calibration of optical detectors for the Yakutsk 
array and for Fly's Eye may result in a systematic disag- 
reement of their data even in the range of high B 0 where 
Fly's Eye geometric factor has been known quite accurately. 
For instance, a Z% variation of the absolute Value of ene- 
rgy leads to a complete agreement between the Yakutsk^yland 
Fly's Eye data at E Q > lO^eV, although it results in an 
increased disagreement between their spectra at E p CiO^eV, 
A more essential disagreement is observed between, on the 
one hand, the data obtained at Haverah Park and the Univer- 
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verelty of Sydney [ll} where no direct energy calibration 
wee mad* and, on the other hand, the data obtained with the 
Yakutsk array and the Fly’s Eye which were energy-callbra- 

ui m K 

The Haverah Park and Sydney arrays give a flat spec- 
trum in the range from E > lO^eV up to the highest detec- 
table energies above 10^ eV * it is not clear whet is the 
role played by the fluctuations of the charged particle la- 
{MMfc W&y distribution function in determining the value of 
$ 600 for ^e Haverah Perk array because the value of 
at E^ > 10^®eV is not determined by interpolation, but is 
found by specifying a fixed mean lateral distribution fun- 
ction. This circumstance may prove to be oven more impor- 
tant for the Sydney array because of a great distance bet- 
ween its detectors. With a foiling ^600 s P ec t**um tor the 
spectrum of , the total number of EAS muons, in case of 
the Sydney array), the great errors In determining the BAS 
sise may flatten the measured spectrum. An additional so- 
urce of errors arises also when we treat the showers whose 
axes fall at t e array periphery. 

So, we think it necessary (1) to rigorously restrict, 
making allowance for the lateral distribution function flu- 
ctuations, the effec tive detec lion area for FAS with E > 

4Q 0 

> 10 eV by the region where the accuracy in determining 
the PAS parameters will make sure that errors would be ab- 
sent and (2) to carry out the Monte Carlo simulation of the 
entire process of measuring and analyzing EAS with an expe- 
rimental array in order to obtain a distsibfton function of 
the type shown in Fig, 5. 

•jkV 11 : """ '* ************* i“m»irirr,,mr,i-,Tf „ rn ■.«!,** — MMMM* 

'The energy spectrum inferred from the Sydney array data 
ie not shown in Pigs 8 and 9 because it is model-dependa 
dent. It is necessary to use the results of the Yakutsk 
array calibration [g] and to scale the muon number spec- 
trum to the primary spectrum. 
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Thus, oh examining the available Yakutsk and Ply's Bye 
data* one may arrive at the conclusion that they do not 
oontradiot qualitatively the Zatsepin-Orelsen pattern of % 
the cutoff of the cosmic ray energy spectrum due to the 
interactions with the 2,71 universal micrawave radiation. 
Additional studies have to be carried out to demonstrate 
the existence of such a cutoff and to find the 'fine* 
spectral structure (bump) , 
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Yakutsk EAS array 
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o - scinti licit con detector f 4m z 
• -scintillation detector, 2 m 2 
v-Cenenfcov tight detector, 
x- Cerenkov pulse width detector, 
500 cm z 

v- Cerenkov pulse width detector, 

ZOO cm 2, 

Q- muon detector, 36 n* 2 
detector, J m 1 
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